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The environments of a complete subsample of 6 of the most powerful radio-loud AGN at 
redshifts z ^ 1.6 are investigated, using deep RJK imaging to depths of R ^ 26, J ^ 22.4 
and K ^ 20.6. An excess of galaxy counts in the K-band is seen across these fields; these 
surplus galaxy counts are predominantly associated with red galaxies {R — K > 4) of magni- 
tudes 17.5 < K < 20.5 found within radial distances of ^1 Mpc of the AGN host. These are 
1^ ' exactly (though not uniquely) the magnitudes, colours and locations that would be expected 

\ of old passive elliptical galaxies in cluster environments at the redshifts of these AGN. Using 

, both an Abell-style classification scheme and investigations of the angular and spatial cross- 

' correlation functions of the galaxies, the average environment of the fields around these AGN 

\ is found to be consistent with Abell cluster richness classes and 1 . The amplitude of the 

CO i angular cross-correlation function around the AGN is shown to be a strong function of galaxy 

' colour, and is highest when only those galaxies with the colours expected of old elliptical 

galaxies at these redshifts are considered. 
pi i' The images cover a relatively wide field, 5 by 5 arcmins, allowing the distribution of the 

surplus galaxy counts to be investigated. The galaxy overdensities are found on two scales 
around the AGN: (i) pronounced central concentrations on radial scales within ^ 150 kpc; 
where present, these are composed almost entirely of red (i? — X > 4) galaxies, suggesting 
^ \ that the morphology-density relation is imprinted into the centres of clusters at a high redshift. 

(ii) weaker large-scale excesses extending out to between 1 and 1.5 Mpc radius. The presence 
or absence of galaxy excesses on these two scales, however, differs greatly between the six 
k> \ different fields: the fields of two AGN do show red galaxy excesses on both scales, another 

j_j ■ two fields show only a large-scale red galaxy overdensity with no pronounced central con- 

' centration, and one field shows only a sharp central peak of red galaxies with no large-scale 

overdensity. The final field shows little evidence for an excess on any scale; this field is as- 
sociated with an unresolved radio source, perhaps indicating that only extended radio sources 
probe cluster environments. 

Clearly there is a large range in both the richness and the degree of concentration of any 
clustering environments around these distant AGN. The implications of this for both cluster 
formation and the nature of high redshift AGN are discussed. 

Key words: galaxies: clusters: general — radio continuum: galaxies — galaxies: evolution 



o 



1 INTRODUCTION 

Understanding the formation of large scale structure and the related 
formation and evolution of galaxies are two of the most important 
issues in modern-day astronomy. Rich clusters of galaxies at high 
redshifts (z > 1) provide an opportunity to investigate both of these 
problems directly: first, since clusters are the largest, most mas- 
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sive, collapsed structures in the Universe, discerning the epoch and 
process of their assembly places strong constraints on models of 
large scale structure formation (e.g. Bahcall & Fan 1998); second, 
because they contain large numbers of galaxies at the same dis- 
tance, clusters provide a unique resource for investigating galaxy 
evolution through the evolution of galaxy scaling relations such as 
the fundamental plane (e.g. van Dokkum et al 1998) and colour- 
magnitude relations (e.g. Stanford et al. 1998); third, clusters con- 
tain the oldest galaxies known, and these can strongly constrain the 
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first epoch of formation of ellipticals (cf Dunlop et al. 1996, Spin- 
rad et al. 1997). 

The identification of significant samples of clusters much 
above redshift one is both challenging and inefficient using opti- 
cal or X-ray selection techniques. Existing X-ray surveys suffer 
from sensitivity limits (the ROS AT Deep Cluster Survey, for exam- 
ple, found very few clusters above z ~ 1, and none above 2 = 1.3; 
e.g. Stanford et al. 2002) whilst the smaller field of view of the 
current generation of X-ray telescopes make them relatively inef- 
ficient for wide-area surveys. At optical wavelengths the contrast 
of a cluster above the background counts is minimal above redshift 
one, requiring very deep wide-area multi-colour surveys including 
near-infrared wavebands to identify a significant sample of 2 1 
clusters; whilst certainly feasible, this is very expensive on tele- 
scope time. For this reason, an alternative method has frequently 
been adopted, using targetted studies towards powerful AGN. At 
low redshifts AGN are frequently found in relatively rich environ- 
ments (e.g. Yates, Miller and Peacock 1989; Hill & Lilly 1991), 
and the host galaxies of high redshift AGN are amongst the most 
massive galaxies known in the early Universe (e.g. de Breuck et al. 
2002), making these promising candidates for residing in galaxy 
overdensities. 

The environments of powerful radio sources at redshifts 2^1 
have been investigated by many authors over recent years using 
a variety of different techniques, including searches for luminous 
extended X-ray emission (e.g. Crawford & Fabian 1996; Fabian 
2001), searches for over-densities of galaxies in infrared colours 
(e.g. Dickinson 1997, Hall & Green 1998) or emission line images 
(e.g. McCarthy et al. 1995), investigations of the colour-magnitude 
relations for red cluster ellipticals (e.g. Best 2000, hereafter BOO), 
cross-correlation analyses (e.g. BOO; Wold et al 2000) and direct 
spectroscopic studies of individual sources (e.g. Dickinson 1997, 
Deltom et al. 1997). The evidence that at least some powerful 2^1 
radio sources are located at the centres of galaxy overdensities is 
overwhelming (see BOO for a review). What remains unclear is the 
ubiquity, scale, and nature of these (proto?) clusters. 

At higher redshifts, radio sources have been detected out to 
z = 5.2 (van Breugel et al. 1999), and some well-studied sources 
have been spectroscopically confirmed to lie in overdense environ- 
ments (e.g. 1138-215 at 2 = 2.2; Kurk et al. 2000; Kurk et al. 
in preparation). An on-going VLT large project (PI: Miley) has 
found order-of-magnitude overdensities of Ly-a emitters around 
all 4 deeply studied 2 > 2 radio galaxies to date (Venemans et al. 
in preparation), including the detection of the most distant known 
large structure around the z = 4.1 radio galaxy TN J1338-1942 
(Venemans et al. 2002). The ~ 20 line emitting galaxies in the 
structures around each of these high redshift radio galaxies are 
found over regions extending to > 3 Mpc, and have velocity dis- 
persions of 300-1000 km s "^(Venemans et al., private communi- 
cation). Thus, it does indeed seem that powerful radio galaxies 
offer unique^ prospects for tracking galaxy clusters over cosmic 
time. The detection of overdensities at these redshifts, however, 
is based predominantly on observing emission-line galaxies near 
the radio galaxy redshift; the strong red sequences characteristic of 
nearby clusters are not observed towards these 2 > 2 overdensi- 
ties, restricting the evolutionary studies of cluster ellipticals using 



t Of course, within a few years deeper X-ray surveys and Sunyaev- 
Zel'dovich effect surveys will also find clusters at redsliifts 2 S> 1, if such 
clusters already contain a significant hot gas component. 



the fundamental plane and colour-magnitude relations from being 
extended to these high redshifts. 

The redshift range 1.2 ^ z ^ 2.0 is therefore a particu- 
larly interesting redshift range to study. This is the redshift range 
where cluster studies go beyond those currently being carried out 
using optical and X-ray selected clusters, whilst still being of low- 
enough redshift that red sequences of cluster ellipticals are likely 
to exist, and can be studied in a practical amotmt of observing 
time. Such clusters will also allow a link to be made between the 
'proto'-clusters being found around 2 > 2 radio galaxies, and the 
rich virialised clusters seen in the nearby Universe. The elliptical 
galaxy population in clusters may well be undergoing a fundamen- 
tal change within this redshift range 

Hall & Green (1998) have previously investigated the small- 
scale environments of a sample of 31 quasars with redshifts 1 < 
2 < 2, variably selected from previous surveys. They found an ex- 
cess of K-band galaxy counts on two scales around the quasars; 
a pronounced peak in the galaxy counts within the central 40 arc- 
sec radius surrounding the quasars was accompanied by a uniform 
galaxy overdensity across their entire fields (to about ~ 1.5 ar- 
cmins radius). Comparing R and K-band data they showed that the 
excess counts correspond to redder galaxies than the typical field 
population, with a significant overdensity of extremely red objects. 
This is indicative that a red elliptical sequence may be present in 
these fields. Hall et al. (2001) investigated the physical extent of 
the large-scale overdensity in the fields of two of these AGN, and 
found a net excess of galaxy counts out to 140 arcsec (~ 1.3 Mpc) 
in one of the two fields. 

We have embarked upon a project to investigate in detail the 
environments of a complete sample of the most luminous radio 
sources (radio galaxies and quasars) within the narrow redshift 
range 1.44 < 2 < 1.7. The aim of this project is to use deep, 
wide-area multi-colour imaging to study in detail any galaxy over- 
density, in particular the nature of the excess galaxies, the presence 
of a red sequence of cluster ellipticals, and the radial extent of the 
overdensity. By studying a complete AGN sample it will be pos- 
sible to investigate the variety of environments found, and whether 
the environment of a radio-loud AGN is in any way correlated with 
the radio source properties. The candidate cluster galaxies in the 
richest overdensities can then be foUowed-up using multi-object 
spectroscopy in order to extend the galaxy evolution studies. 

This paper provides the initial results of this project, focussing 
on the imaging data and results from the first 6 clusters in the sam- 
ple. Specifically, the lay-out of the paper is as follows. In Sec- 
tion 2 the sample selection, observations, and data reduction are 
described. The resulting multi-colour images of each field are pro- 
vided and described in Section 3, and galaxy catalogues are con- 
structed from these. The galaxy number counts and the colour dis- 
tribution of the galaxies are investigated in Section 4, and a cross- 
correlation analysis carried out. The properties of the different 
fields are related to the properties of their AGN in Section 5. The 
results are discussed and conclusions drawn in Section 6. Through- 
out the paper, the values adopted for the cosmological parameters 
are CIm = 0.3, f^A = 0.7 and Ho = 65kms"^Mpc"\ 



2 OBSERVATIONS 

2.1 Sample Selection and Observations 

The radio sources were drawn from the equatorial sample of 178 
powerful radio sources defined by Best, Rottgering and Lehnert 
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(1999; the BRL sample). This sample was designed to be equiv- 
alent to the northern 3CR sample, and was drawn from the Molon- 
glo radio survey at 408 MHz (Large et al. 1981), according to three 
criteria: (i) the sources must have declinations —30° < 5 < +10°; 
(ii) they must lie away from the galactic plane, |6| > 10; (iii) they 
must be brighter than 5'408mhz = 5 Jy (roughly equivalent to the 
definition of the 3CR sample). Spectroscopic redshifts are available 
for 177 of these 178 sources (Best et al. 1999; Best et al. 2000). 

From the BRL sample, the 9 radio sources with redshifts 
1.44 < z < 1.7 and galactic latitudes |6| > 25 were selected. 
This sample comprised 4 radio galaxies and 5 radio loud quasars. 
Of these, the environments of 6 radio sources have been studied to 
date, and are presented in this paper. These 6 radio sources, limited 
by the constraints of telescope time, correspond to all of those with 
15 < RA < 02 and therefore constitute a complete subsample of 
the most powerful radio sources at these redshifts. 

The aim of the current observations was to deeply image these 
fields in the R, J and K wavebands. This filter combination was 
chosen to provide maximal sensitivity to old cluster ellipticals: the 
strongest continuum features of such ellipticals is the 4000A break, 
which at z ~ 1.5 falls between the R and J bands, making this 
colour very sensitive to old cluster galaxies. The addition of the 
Ti'-band data provides a long red colour baseline to sample the 
old stellar populations of these ellipticals, whilst also facilitating 
identification of bluer cluster galaxies with on-going star formation 
(and hence less pronounced 4000A breaks), and helping to filter 
out non-cluster objects with very red R — J colours, such as dusty 
starbursting galaxies. At z ~ 1.5 passively evolving old ellipticals 
would have colours of J — if ~ 1.8 and R—J^3.7, and brightest 
cluster galaxies at these redshifts are expected to have magnitudes 
of K ^ 17.5 to 18 (cf. those of radio galaxies at these redshifts, 
e.g. Best et al. 1998). The goal of detecting elliptical galaxies in all 
colours to ~ 3 magnitudes fainter than this set the required depths 
of the different observations. Full details of the observations are 
provided in Table 1. 

The six fields were observed in the J and if -bands using 
the infrared imager SOFI on the ESO New Technology Telescope 
(NTT) during September 2000 and August 2001. SOFI has a 1024 
by 1024 HgCdTe array with 0.29 arcsec pixels, providing a field of 
view of 5 by 5 arcmins. Images were taken using a jittering tech- 
nique, with the pointing centre moved randomly within a 40 arcsec 
square box between each exposure. At each pointing position, 1- 
minute of data was taken, split into six 10s co-adds in the K-band 
and three 20s co-adds in the J-band. In the if -band the Kg filter 
was used to reduce the sky background; the Kg magnitude is re- 
lated to the if -magnitude by K - Ks ~ -0.005(J - if), and 
all quoted magnitudes have been converted to the if -band mag- 
nitude. Between 1.7 and 2.5 hours of data were obtained for each 
source (see Table 1 for details), reaching a 'limiting magnitude' of 
if ~ 20.6; this value corresponds to the 50% completeness limit 
for galaxies with sizes typical of those at redshifts z ~ 1.6. In 
the J-band the total exposure time varied between 1 and 2 hours, 
reaching a limiting depth of J ~ 22.4. 

For 5 of the 6 fields (2128—208 was randomly excluded due 
to telescope time constraints), i?-band data for the field were taken 
through the Bessel R filter using SUSI-2 on the NTT. SUSI-2 com- 
prises two 2048 by 4096 pixel EEV CCDs, with a pixel scale of 
0.08 arcsec, which combined provide a 5.5 by 5.5 arcminute field 
of view. The two CCDs are separated by about 10 arcseconds on 
the sky; the total exposure time was therefore separated into 5-min 
blocks, with the pointing position moved through a 9-point pat- 
tern of offsets [(0,0), (10,10), (20,-10), (30,0), (40,10), (-40,-10), (- 



30,0), (-20,10), (-10,-10) arcsec] in order to attain equal sensitivity 
over the entire SOFI field. This pattern was repeated until between 
3.75 and 6 hours of data had been obtained on each field, reaching 
a limiting depth of i? ~ 25.9. 

Standard stars were observed approximately every 2 hours 
throughout the nights to allow the images to be photometrically 
calibrated. For the SUSI-2 images, once each night the image was 
rotated by 1 80 degrees in order to place the standard stars on the op- 
posite CCD; no significant differences in calibration were apparent 
between the two CCDs. All of the nights were photometric, and the 
instrument zero-points varied by less than 0.05 magnitudes within 
each of the Sept 2000 and Aug 2001 runs, although the i?-band 
zero point changed by 0.1 magnitudes between the two runs. 

2.2 Data reduction 

The infrared observations were reduced in the standard manner, us- 
ing the DiMSUM^ package within iraf. After subtraction of a dark 
frame, and removal of known bad pixels, a median of the images 
(after masking of bright objects) was constructed to obtain an ac- 
curate sky flat-field. The flat-fielded images were then sky sub- 
tracted, cosmic ray events were removed, and then the individual 
frames were block replicated by a factor of 2 in each dimension to 
allow more accurate alignment, resulting in a final image scale of 
0. 144 arcsec per pixel. The frames were accurately registered using 
the peak positions of several bright unresolved objects visible on all 
images, and were combined. Because of the long exposure times, 
the seeing sometimes varied throughout the dataset, and so during 
the combination process the images were weighted according to 
their seeing in order to provide a final image with the optimal com- 
bination of seeing and sensitivity. Stars from the US Naval Obser- 
vatory catalogue A2.0 (Monet et al. 1998) were used to provide an 
astrometric solution for the images. The absolute astrometry should 
be accurate to better than 0.5 arcsec. The fields were photometri- 
cally calibrated onto the standard Vega photometric system. 

The ii-band images were reduced within IRAF. The images 
were bias subtracted, and then flat-fielded using a flat-field con- 
structed by from the median of the various images (after masking 
bright objects). An illumination correction was then applied to re- 
move large-scale variations in the flat-fielded images. The images 
were registered using the peak positions of several bright unre- 
solved objects, and then they were combined, weighted according 
to their seeing as for the J and if band data. Cosmic rays were 
rejected during the combination process. The resultant image was 
then astrometrically transformed using a 5 parameter fit (RA and 
Dec of the central pixel, pixel scale in x and y directions, and rota- 
tion angle) to align it with the infrared data, converting it to 0. 144 
arcsec pixels whilst conserving the object magnitudes. The relative 
astrometry of the IR and optical frames is accurate to < 0.1 arcsec. 



3 THE RADIO SOURCE FIELDS 

K-band images of the six fields are shown in Figures 1 to 3, with 
the radio source host galaxy marked by a cross in each case. Three- 
colour images of the 100 by 100 arcsec region around the radio 
source host in each field are provided in Figures 4 to 6. A brief 
discussion of each individual field follows. 

t DiMSUM is the 'Deep Infrared Mosaicing Software' package developed 
by Eisenhardt, Dickinson, Stanford and Ward. 
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Table 1. Details of the observations. 

Source RA Dec z Waveband Obs. Telescope Exp. Mean Limiting 

J2000 Dates & Time Seeing Magnitude* 

Instrument [mins] [arcsec] 



MRC 0000-177 00 03 21.93 -17 27 11.9 1.47 

J 
R 

MRC 0016-129 0018 51.36 -UAliAA 1.59 K, 

J 
R 

MRC 0139-273 01 41 27.31 -27 06 11.1 1.44 K, 

J 
R 

MRC 1524-136 15 26 59.45 -13 5100.2 1.69 K, 

J 
R 

MRC 2025-155 20 28 07.68 -15 21 20.9 1.50 K, 

J 
R 

MRC 2128-208 21 31 01.49 -20 36 56.4 1.62 K, 

J 



2001 Aug 10,12 


NTT-SOFI 


138 


0.98 


20.4 


2001 Aug 11,12 


NTT-SOn 


120 


1.01 


22.3 


2001 Aug 12 


NTT-SUSI2 


270 


1.00 


25.8 


2000 Sep 23 


NTT-SOFI 


103 


0.62 


20.7 


7000 Spn 1^ 9^ 


NTT-SOFT 


60 


80 


22.3 


2000 Sep 20,24,25 


NTT-SUSI2 


225 


0.85 


25.8 


2000 Sep 23,25 


NTT-SOFI 


135 


0.64 


20.6 


2000 Sep 23,25 


NTT-SOFI 


105 


0.74 


22.6 


2000 Sep 20,24,25 


NTT-SUSI2 


360 


0.74 


26.1 


2001 Aug 10,11 


NTT-SOFI 


135 


0.80 


20.5 


2001 Aug 10 


NTT-SOFI 


105 


0.90 


22.4 


2001 Aug 11,12 


NTT-SUSI2 


255 


1.15 


25.5 


2000 Sep 20,23 


NTT-SOFI 


114 


0.75 


20.4 


2000 Sep 23 


NTT-SOFI 


91 


0.90 


22.3 


2000 Sep 24,25 


NTT-SUSI2 


260 


0.65 


25.9 


2001 Aug 10 


NTT-SOFI 


150 


0.77 


20.6 


2001 Aug 10 


NTT-SOFI 


105 


1.00 


22.5 



* The 'limiting magnitude' quoted is the approximately 50% completeness Unlit for galaxies with physical sizes typical of those ma z ^ 1.6 cluster 



0000-177 (z = 1.47): This radio-loud quasar lies in a relatively 
sparse region of sky. A faint object is seen 1-2 arcsec to the north- 
east of the quasar host, and a further seven objects with R — K>4 

lie within 30 arcsec (~ 275 kpc at z ^ 1.5), indicating a possi- 
ble compact group of galaxies at the quasar redshift. There is no 
obvious large-scale galaxy overdensity. 

0016-129 (z = 1.59): This radio galaxy is the westernmost of a 
string of three luminous red galaxies in the centre of this field, and 
shows extended blue emission along the radio axis direction, char- 
acteristic of the alignment effect in high redshift radio galaxies (cf. 
McCarthy et al. 1987, Chambers et al. 1987). The redder galaxies 
in the field are more concentrated towards the central regions of the 
field, within 1-2 arcmins of the AGN. 

0139-273 (z = 1.44): This is the richest field of those imaged in 
the current observations, and contains a large number of very red 
galaxies. These galaxies are found predominantly within the inner 
regions of the field, although no red galaxies are seen close to the 
AGN. Very luminous elongated blue emission extends from the ra- 
dio galaxy host, along the direction of the radio source axis. 
1524-136 (z = 1.69): This radio loud quasar is the highest red- 
shift AGN in the current sample, and its surrounding field is rel- 
atively rich in the K-band. A moderate excess of red galaxies is 
seen on scales out to about 120 arcsec (~ 1 Mpc). No sharp central 
peak (such as that seen for 2025—155) is seen, suggesting that any 
cluster enviroimient may still be forming. 

2025-155 (z = 1.50): A very compact group of several red galax- 
ies is found within 15-20 arcsec (150 kpc) radius of this radio-loud 
quasar. All of these objects have colours consistent with old passive 
ellipticals at the radio source redshift. A significant overdensity of 
red galaxies is found extending over the entire 5 arcmin field, al- 
though predominantly within the inner regions. 
2128-208 (z = 1.62): This field does not contain a particular ex- 
cess of K-band galaxies and, except for a companion object to the 
quasar, no other galaxy is found within 15 arcsec of the AGN host. 
With only J and K-band data available, there is little information 
about the colour distribution of the galaxies across the field. 

It is clear from this direct examination of the different images 
that a wide variety of environments in fotmd amongst the sample. 



At one extreme, the field of 2128—208 shows little evidence for 
any net overdensity of red galaxies on either small (< 150 kpc) 
or large (~Mpc) scales, whilst in contrast 2025—155 has a very 
rich environment on small scales coupled with an apparent excess 
of red galaxies across the entire field. Between these extremes, the 
other fields show evidence for structure on differing scales: a small 
scale excess of red galaxies is found around 0000— 177, but not 
much large scale structure; the 0016—129 field has three central 
red galaxies and a weak larger-scale excess; above average num- 
bers of red galaxiess are seen on both moderate and large scales 
around 1524—136; and the field of 0139—273 shows a strong ex- 
cess of red galaxies on larger scales, but with few red objects close 
to the AGN. In the discussions of the statistical properties of the 
fields that follow, it is important to bear in mind that global results 
quoted correspond to the average values across 6 rather different 
environments. 

3.1 The source catalogues 

Source extraction and catalogue construction followed the method 
described in detail by BOO, and the reader is referred to that paper 
for full details. This method is summarised briefly below. 

Throughout these analyses, the K-band data were used as the 
primary dataset from which the catalogue is defined. Source detec- 
tion and photometry were carried out on the K-band frames using 
SEXTRACTOR version 2.2.1 (Berlin & Arnouts 1996), using the 
output exposure map produced by dimsum as a weight map for 
SEXTRACTOR to Compensate for the varying noise levels associ- 
ated with the jittering procedure. The source extraction parameters 
were set such that, to be detected, an object must have a flux in 
excess of 1.5 times the local background noise level over at least 
N connected pixels, where N was varied a little according to the 
seeing conditions, but in general was about 15 (equivalent to 4 pix- 
els prior to the block replication during the mosaicing procedure). 
A search for negative holes using the same extraction parameters 
(cf BOO) found only one negative feature over all 6 frames brighter 
than the 50% completeness limit, implying that essentially all of 
the extracted objects down to that level are real. The sextractor 
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Figure 1. [See attached jpg files] Top: The K-band image of the field of 0000-177 (z = 1.47). Bottom: The K-band image of the field of 0016-129 
(z = 1.59). 



Figure 2. [See attached jpg files] Top: The K-band image of the field of 0139-273 (z = 1.44). Bottom: The K-band image of the field of 1524-136 
(z = 1.69). 



catalogues were examined carefully, and obvious problems, such as 
objects associated with the diffraction spikes of bright stars, were 
corrected. Those objects within the central 2000 by 2000 pixels 
(290 by 290 arcsec) of each image were retained for the subse- 
quent analysis; more distant objects lie in the regions of the image 
where the noise levels are higher due to the jittering process, and 
were discarded; at the edges of the retained regions of the fields, 
the noise level is only ~ 10% higher than in the field centre. Small 
regions of sky (of 10, 20 and 10 arcsecond radius) were also dis- 
carded around the three bright stars in the images of 0016—129, 
1524—136 and 2025—155. The six fields combined provide a sky 
area of 140 square arcminutes. 

SEXTRACTOR's MAG_BEST estimator was used to determine 
the magnitudes of the extracted sources; this yields an estimate 
for the 'total' magnitude using Kron's (1980) first-moment algo- 
rithm, except if there is a nearby companion which may bias the 
total magnitude estimate by more than 10% in which case a cor- 
rected isophotal magnitude is used instead, sextractor isopho- 
tal apertures smaller than 1 arcsecond were replaced by circular 
1 arcsecond apertures. The determined magnitudes were corrected 
for galactic extinction using the E{B — V) values provided in 
the NASA/IPAC Extragalactic Database (NED). The accuracy of 
these total magnitudes and the completeness level of the source 
extraction were investigated as a function of position on each im- 
age through a Monte-Carlo simulation, as fully described by BOO. 
Briefly, a total of 10000 stellar objects and 25000 galaxies (see BOO 
for the prescription by which these were constructed) were added, 
25 at a time, to each image with a magnitude between 15 and 22.5, 
and were then re-extracted using sextractor with the same in- 
put parameters as for the original source extraction, to see if the 
added objects were detected. For each recovered source, the differ- 
ence between the total magnitude measured by sextractor and 
the input total magnitude was determined. 

From these results, the mean coinplctcncss fraction was de- 
termined as a function of magnitude for both stars and galaxies 
over the entire set of images. The mean 50% completeness limits 
are K = 20.55 for galaxies and K = 21.05 for stars (see Fig- 
ure 7); the galactic completeness limits for each individual field are 
provided in Table 1. The mean difference between the input and 
measured magnitudes, and the scatter of the measured magnitudes 
aroimd this mean, are shown as functions of input magnitude for 
both model stars and galaxies in Figure 8. This difference was used 
to correct the observed magnitudes and obtain true total magni- 
tudes. 

The J and R band magnitudes for the catalogue sources were 
then measured by running SEXTRACTOR in its double image mode, 
where one image (the K band) is used to detect the objects and the 
fluxes and magnitudes are then measured from the second (J or 
R band) image through the same aperture. Those objects whose J 
or R band fluxes were measured to have a significance of less than 
3a, where a is the uncertainty on the flux measurement provided by 
SEXTRACTOR, were replaced by 3a upper limits. The sextrac- 
tor flux error estimate includes the uncertainty due to the Poisson 
nature of the detected counts and that from the standard deviation 
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Figure 7. The completeness fraction for both stellar (dashed line) and galac- 
tic (solid line) model objects as a function of input magnitude, averaged 
over all 6 fields. 
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Figure 8. The mean (upper plot) and rms scatter (lower plot) of the differ- 
ence between the measured magnitudes and input magnitudes for the model 
stars (dashed line) and galaxies (soUd line). 
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Figure 3. [See attached jpg files] Top: The K-band image of the field of 2025- 
(z = 1.62). 



155 {z = 1.50). Bottom: The K-band image of the field of 2128-208 



Figure 4. [See attached jpg files] Top: An R,J,K colour image of the central 1 00 X 1 00 arcsec field around the z - 
colour image of the central 100 x 100 arcsec field around the z = 1.59 radio galaxy 0016—129. 



1.47 quasar 0000-177. Bottom: An R,J,K 



of the background counts. An additional source of flux error arises 
from the uncertainty in the subtraction of the baclsground count 
level as a function of position across the image. This value was es- 
timated as the product of the area of the extraction aperture and the 
rms variation of the subtracted background flux across blank re- 
gions of the image. This background subtraction error estimate was 
combined in quadrature with the flux error given by SEXTRACTOR 
to determine the uncertainties on the magnitudes of the extracted 
objects in all three colours. 



3.2 Star-Galaxy separation 

SEXTRACTOR provides a 'stellaricity index' for each object, where 
in the ideal case a galaxy has a stellaricity index of 0.0 and a star 
has 1.0; low signal-to-noise, and galaxies more compact than the 
seeing, lead to an overlap in the calculated stellaricity indices for 
the two types of object at the faintest magnitudes, and hence uncer- 
tain star-galaxy separation. Stellaricity indices were taken from the 
K-band data. At magnitudes K < 18, star-galaxy separation could 
be carried out directly from these stellaricity values taking values 
above 0.8 (the generally adopted cut-off) to correspond to stars. At 
fainter magnitudes, separation is less certain, but was aided by the 
addition of the multi-colour information; by plotting the R-K ver- 
sus J-K colour-colour diagram as a function of stellaricity index, a 
region of colour space typically occupied by stars was defined (see 
Figure 9), and all objects in this colour space region with stellar- 
icities above 0.5 were also classified as stars. This classification is 
somewhat conservative, and a small proportion of stars may be in- 
cluded along with the galaxies, but this was considered preferable 
to falsely classifying galaxies as stars. 

Figure 10 shows the total star counts as a function of mag- 
nitude derived from all of the frames, fit with the function 
log N{K) = 0.201/s: - 0.128. This relatively good fit to the dif- 
ferential star counts using a single power-law distribution demon- 
strates that the star-galaxy separation is working well to at least 
K ~ 20. All objects classified as stars were removed from further 
analysis. 



4 GALAXY COUNTS AND COLOURS 
4.1 Galaxy counts 

The galaxy counts were determined by once again following the 
procedure of BOO. In brief, the measured magnitudes of the galax- 
ies were converted to true 'total' magnitudes using the offsets de- 
termined from the simulations above, and binned in 0.5 magnitude 
bins over the magnitude range 14 < K < 21. These raw galaxy 
counts, riraw, were then corrected for completeness and for magni- 
tude biasing effects (again as calculated from the results of the sim- 
ulations) to produce corrected galaxy counts (ric). Scaling by the 
observed sky area produced final counts per magnitude per square 
degree, Nc. In calculating the uncertainty on this value, in addi- 
tion to the Poissonian error on the raw galaxy counts, an extra error 
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Figure 9. A set of i? — X versus J — K colour-colour plots for objects 
in the field of 2025-155, separated by the stellaricity index output by SEX- 
TRACTOR. The solid lines define a region of colours separating the majority 
of the stellar locus from typical galactic colours. AH objects with stellaricity 
indices greater than 0.8 were classified as stars, together with those objects 
with stellaricities between 0.5 and 0.8 lying to the upper left of the two 
lines. 
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Figure 10. The differential star number counts in the K-band, in AK 
0.5 magnitude intervals, fit using a single power law distribution. 



term was included to account for the uncertainty of the complete- 
ness correction procedure; this was estimated at 20% of the number 
of galaxies added in the correction. The resulting galaxy counts are 
tabulated in Table 2 and are plotted in Figure 1 1 compared with 
counts from various K-band blank-field surveys (Gardner et al. 
1993; Mcleod et al. 1995; Djorgovski et al. 1995; IVIoustakas et al. 
1997; Minezaki et al. 198; Szokoly et al. 1998; Bershady et al. 
1998; Saracco et al. 1999; Vaisiinen et al. 2000; Totani et al. 2001). 
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Figure 5. [See attached jpg files] Top: An R,J,K colour image of the central 100 x 100 arcsec field around the z = 1.44 quasar 0139—273. Bottom: An R,J,K 
colour image of the central 100 x 100 arcsec field around the z = 1.69 radio galaxy 1524—136. 



Figure 6. [See attached jpg files] Top: An R,J,K colour image of the centra! 100 X 100 arcsec field around the z = 1.50 quasar 2025—155. Bottom: A J,K 
colour image of the central 100 X 100 arcsec field around the z = 1.62 quasar 2128—208. 
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Figure 11. The K-band galaxy number counts per magnitude per square 
degree. The solid symbols represent the data presented in this paper and 
tabulated in Table 2. The open symbols represent data from the litera- 
ture: squares are from Minezaki et al. (1998), triangles from Djorgovski 
et al. (1995), diamonds from Moustakas et al. (1997), circles from Szokoly 
et al. (1998), x's from Gardner, Cowie & Wainscoat (1993), asterisks from 
McLeod et al. (1995), inverted triangles from Bershady, Lowenthal & Koo 
(1998), stars from Totani et al. (2001), crosses from Vaisanen et al. (2000), 
and bow-ties (M) from Saracco et al. (1999). 

The derived counts are fully consistent with the literature 
counts at magnitudes K < 17.5. This magnitude corresponds to 
the expected magnitude of brightest cluster galaxies at z ~ 1.5. 
Fainter than this magnitude an surplus of galaxies is seen; the 
galaxy counts still lie within the scatter of the other observa- 
tions, that scatter being largely dominated by cosmic variance (cf 
Daddi et al. 2000a), but the counts in every magnitude bin from 
17.5 < K < 20.5 exceed the literature average by at least la. 

In Table 3 the total galaxy counts in the magnitude range 
17.5 < K < 20.5 are provided for each individual field, and com- 
pared to the mean value of the blank-field surveys. To carry out 
such a comparison it is important to bear in mind that the uncer- 
tainty on the counts in the AGN fields will be significantly larger 
than solely the Poissonian error, due to cosmic variance and the ef- 
fects of galaxy clustering. The true tmcertainty on the counts can 
be estimated from the equation (J^,-uc = -^(1 + NAC), where A'^ 
is the total counts, A is the angular cross-correlation amplitude, 
and C is the integral constraint which depends upon the size and 
geometry of the field. As derived in Section 4.6, for these fields 
the values of A and C are 0.0005 and 16.7 respectively. Hence, for 
N ~ 300 galaxies detected in each field, the uncertainty on the 
counts is approximately double that which would be expected from 
Poissonian statistics alone. The errors tabulated in Table 3 are these 
true errors, including this contribution of cosmic variance. For the 
literature counts, given that these are derived by combining many 
independent surveys over different areas of slsy, the combination 
should average out the cosmic variance effects to a large extent. 

It is noteworthy that five of the six fields have counts above 
the literature average, with three of the fields being > 2cr above. 
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Table 2. K-band galaxy counts as a function of magnitude, over the com- 
bined 140 square arcminutes. The columns give the raw (riraw) and cor- 
rected (ric) counts, the corrected counts per square degree per unit magni- 
tude (Nc), the error on this, and the mean literature coimts (Nm). 
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20.0 
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23100 
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3500 
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21.0 


292 


927.8 


47900 
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41800 



* : counts per magnitude per square degree. 

Table 3. Properties of the galaxies in the individual fields. The columns 
give: (1) the radio source field; (2) the K-band galaxy counts (in galaxies 
per square degree) between 17.5 < K < 20.5, and the uncertainty of this 
value — this uncertainty includes both the Poissonian error and a cosmic 
variance term; (3) the deviation of the number counts from the mean lit- 
erature counts; (4 and 5) the mean and median R — K colours of these 
galaxies; (6) the AbeU cluster strength estimator iVo.s, and its error. The 
bottom two rows of the table give the values average over all of the fields, 
and blank-field expectations for the number coimts and average colovirs as 



drawn from the Uterature (see text for details) 
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[deg-2] 




Mean 


Median 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0000-177 


42200 ± 5100 


-l.Oo- 


4.19 


4.19 


-3 ±5 


0016-129 


55400 ± 5800 


+ lAa 


4.09 


4.05 


4±6 


0139-273 


64200 ± 6200 


+2.7cr 


3.98 


3.93 


18 ± 8 


1524-136 


59600 ± 6100 


+2.0(7 


4.07 


3.85 


14± 7 


2025-155 


58700 ± 5900 


+ 1.9cr 


4.05 


3.89 


11 ± 7 


2128-208 


50400 ± 5600 


+0.6cr 






3±6 


Combined 


55000 ± 2400 


+3.1<T 


4.08 


3.95 


8±3 


Literature 


47500 ± 2000 




3.87 


3.78 





Combining the fields, the cumulative counts over this magnitude 
range are 3a above the average of blank-field surveys; it should be 
emphasised that these quoted significance levels including the ef- 
fects of cosmic variance, and are twice this size considering Pois- 
son statistics alone. Indeed, since the cosmic variance is caused by 
the clustering of galaxies, fields which contain clusters of galaxies, 
like these are proposed to, are in fact the type of field that con- 
tributes significantly to this cosmic variance, and so the Poissonian 
estimate may be more appropriate. 

The cumulative excess of galaxy coimts corresponds to about 
290 (±100) additional galaxies, that is, of order 50 galaxies per 
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field. Although it is possible that some fraction of these correspond 
to faint stars that were misclassified as galaxies (cf. the fall in star 
counts in the final bin of Figure 10), it is clear that these fields are, 
on average, richer than average environments. For comparison, in 
Abell's original cluster definition (Abell 1958) he considered the 
excess counts of galaxies within a radius of \.5h~^Mpc of the 
cluster centre, with magnitudes between ms and ms + 2, where 
mz is the magnitude of the third ranked cluster member: an ex- 
cess of 30 to 49 galaxies corresponded to Abell Class clusters, 
and an excess of 50 to 79 galaxies was Abell Class 1. Assuming 
that the radio host galaxies are of order a magnitude brighter than 
a third-ranked cluster member§ the average excess number counts 
correspond to moderate richness clusters. 

To obtain a better indication of the cluster richnesses, the pa- 
rameter ATo.s defined by Hill and Lilly (1991) can be used. This 
is an Abell-type measurement defined as the net excess number 
of galaxies within a projected radius of 0.5 Mpc (at the cluster 
redshift) of the central galaxy with magnitudes between mi and 
mi -|- 3, mi being the magnitude of the brightest cluster galaxy^. 
The value of iVo.s was calculated for each field, in each case using 
the region of the image at projected distances greater than 1 Mpc 
distant from the AGN to estimate the background counts, thus ac- 
counting to a large an extent as possible for the cosmic variance in 
the background sources. 

The values of A^o.s derived are provided in Table 3. Although 
the errors on each individual value are large, it is notable that 
the three values significantly above zero correspond to the three 
fields with large excesses in their number counts. This implies that 
the number count excesses are predominantly in the iimer regions 
around the AGN, and that these are therefore consistent with clus- 
ter environments. Hill and Lilly (and see also Wold et al. 2000) 
calculated a conversion between the value of Afo.5 and the tradi- 
tional measure of Abell richness. N0.5 ~ 10 corresponds to the an 
Abell richness cluster, whilst A^o.r, ~ 20 corresponds to Abell 
richness 1. Three of the 6 fields studied are compatible with being 
relatively rich clusters, whilst the other three seem to be at most of 
Abell class richness. It should be noted, however, that conclusions 
based upon galaxy counts alone are not robust. 

4.2 Galaxy colours 

The presence of red galaxies associated with a cluster ought to red- 
den the average galaxy colours with respect to blank-field expec- 
tations, and indeed this is seen in the mean or median colours of 
each field (see Table 3). Interestingly, the average colour in the 
radio galaxy fields is typically reddest in the fields with the low- 
est number counts, whereas naively the opposite result would be 
expected if the excess number counts were solely associated with 
typically red cluster galaxies. This result may imply that there are 
excess red galaxies around all of the AGN, and in the cases where 
the non-cluster (foreground and background) galaxy counts are ran- 
domly lower, the red cluster galaxies are therefore a larger propor- 

§ If the radio galaxies were dramatically brighter in the K-band than other 
cluster members, as is the case for some cD galaxies in nearby clusters, then 
that would require the fields to be very rich, probably implausibly so. 
^ For the quasars, the value of rn\ was estimated from the redshift using 
the Hubble K— 2; relation for radio galaxies of the same radio power (e.g. 
LUly & Longair 1984, Best et al. 1998); this is appropriate due to the simi- 
larity of radio galaxy and quasar hosts at high redshifts (e.g. Lehnert et al. 
1992, Kukulaetal. 2001). 



Table 4. A comparison of the median R — K colours (together with the 
errors on the mean) of galaxies in the central regions of the AGN fields, the 
entirety of the AGN fields, and the blank-field distribution from the CADIS 
16hr field and K20 survey. 
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16-17 


3.44 


0.13 


3.50 


0.06 


3.45 


0.07 


17-18 


3.69 


0.14 


3.75 


0.07 


3.65 


0.05 


18-19 


4.12 


0.11 


4.02 


0.06 


3.83 


0.04 


19-20 


4.20 


0.09 


4.08 


0.05 


3.86 


0.03 


20-21 


3.98 


0.06 


3.89 


0.03 


3.75 


0.02 



tion of the total counts and influence the average colours more. This 
should offer further caution about interpreting the cluster environ- 
ments based on counts in a single colour only. 

To investigate in more detail the nature of the excess galax- 
ies, the R— K colour distribution of the galaxies was derived as a 
function of projected distance from the AGN for all galaxies with 
K-magnitudes in the range 16 < /sT < 21 across the five fields 
(excluding 2128—208 for which J?-band data is lacking). Four 
radial bins were considered (0-100, 100-250, 250-500 and 500- 
1000 kpc, as evaluated at the redshifts of the AGN), and the colour 
distribution of the galaxies in these radial ranges are shown in Fig- 
ure 12. Also shown as a dotted line in these plots is the colour 
distribution of those sources beyond 1 Mpc radius from the AGN, 
constructed as a comparison sample. There is a clear surplus of red 
(i? — if > 4) galaxies at all radii < 1 Mpc compared to the outer 
'blank-field' regions. This is most pronounced in the inner 100 kpc 
radius, but is very significant all the way out to 1 Mpc. 

This colour distribution was also derived as a function of K- 
band magnitude, both over the total combined area of the five im- 
ages and over just the inner quadrant of the images, correspond- 
ing to distances out to about 700 kpc from the AGN. These dis- 
tributions are displayed in Figure 13. The median colours of the 
galaxy distributions in each magnitude range can be found in Ta- 
ble 4. To compare these, a blank-field distribution was also deter- 
mined by combining the data of the CADIS 16hr field (Thomp- 
son et al. 1999; data kindly supplied by Dave Thompson) and the 
K20 survey (Cimatti et al. 2002; data kindly supplied by Andrea 
Cimatti). The data from the CADIS survey were colour corrected 
to account for differences in the R and K passbands according 
to the relation R - K = {R - K')cAms + 0.07. The distribu- 
tions of the CADIS and K20 data were then well matched. These 
two fields cover a combined sky area of 206 square arcminutes, 
and are complete for K < 20. For 20 < < 21, a blank-field 
galaxy R — K distribution was constructed using the compilation 
of Hall & Green (1998), colour correcting the passbands according 
\.oR- K ={r- K)iic, - 0.26 (Hall & Green 1998). These ex- 
pected blank-field distributions are indicated by the shaded regions 
in Figure 13 and the median colours tabulated in Table 4. 

At magnitudes K < 18, the colour distribution of the total 
fields and the central regions are both well matched to the field dis- 
tribution, with comparable median colours. By 18 < if < 19 the 
median colour begins to become redder in the AGN fields than in 
the blank fields, particularly in the central regions, due to a small 
excess of galaxies with R — K > 5. This red excess becomes even 
more pronounced in the two fainter bins. The combined colour dis- 
tribution for all galaxies with 16 < if < 21 clearly shows that the 
higher number counts in these fields are associated with red galax- 
ies: the blank-field predictions are matched very well for galaxies 
bluer than J? — if = 3, but the AGN fields show a strong surplus 
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Figure 12. The R— K colour distribution of tlie galaxies as a function of radius, combined over the 5 fields. The solid lines show the colour distribution in 
each range of projected distances from the AGN (as evaluated at each radio source redshift), and the dashed lines show the distribution constructed from the 
regions of the images at projected radial distances greater than 1 Mpc. A large excess of red galaxies is seen, most pronounced in the inner 100 kpc but also 
very significant even out to distances of a Mpc. 



of redder galaxies, particularly those with R — K > 5. These re- 
sults are even more pronounced when only the central regions of 
the images, around the AGN, are considered. 

Combining the results for the colour distributions of the galax- 
ies as a function of both radius and magnitude, it is clear that the 
predominant excess galaxy population corresponds to a population 
of red (R — K > 4 — 5) galaxies, with magnitudes K > 18, lying 
at projected distances within ~ 1 Mpc from the AGN. These are 
exactly the typical magnitudes, colours and locations that would be 
expected of elUptical galaxies in a cluster environment at the AGN 
redshifts, if such ellipticals form at early cosmic epoch and evolve 
relatively passively. By contrast, in semi-analytic models of galaxy 
formation, very few ellipticals this red and this luminous would be 
expected in the field at redshifts a ~ 1.6 (cf Daddi et al. 2000b). 

4.3 Radial distribution of individual fields 

Knowing the magnitudes and colours of the galaxies producing the 
excess counts, it is possible to investigate their radial distribution in 
each of the galaxy fields separately. Considering only those galax- 
ies with magnitudes 17.5 < K < 20.7 and colours R — K > A, 
the radial distribution of the galaxy counts was derived for each 
field. These results are shown in Figure 14, together with a com- 
bined histogram produced by averaging across all five fields. Also 
included on these plots is the expected number density for blank- 
field sources of these magnitudes and colours; this was calculated 



by taking the K20 and CADIS surveys, making the same cuts in 
colour and magnitude as for the AGN fields, and then convolving 
the results with the completeness versus magnitude function (Fig- 
ure 7). 

The results from combining the five AGN fields are in agree- 
ment with those derived by Hall & Green (1998), namely that a 
pronounced peak of red galaxies is seen within the inner ~ 1 50 kpc, 
together with a weaker large-scale excess. With the larger field- 
of-view of the data provided here, it is now possible to quantify the 
radial extent of this large-scale overdensity: the galaxy counts re- 
main above field expectations out to between 1 and 1.5 Mpc radius. 

Comparing the five individual fields in Figure 14, however, it 
is immediately apparent that this combined result is the average of 
five fields with very different properties. 0016—129 and 2025—155 
do both show radial distributions similar to the average result, with 
pronounced peaks in the inner bin followed by weak red galaxy 
overdensities to ~ IMpc. 0000—177 has a central pronounced 
peak of red galaxies, but shows no large-scale excess (indeed quite 
the opposite); this may be indicative of a compact group of galax- 
ies. 0139—273 and 1524—136 are both much less centrally concen- 
trated, with smoother central peaks and much more gradual radial 
declines in red galaxy surface density out to the blank-field ex- 
pectations at radial distances of between 1 and 1.5 Mpc. There is 
clearly a large range in the nature, scale and significance of any 
galaxy overdensities aroimd these AGN. 

That the radial distribution of galaxies within these large-scale 
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Figure 13. The R — K colour distribution of the galaxies as a function of galaxy K-haod magnitude. The solid lines represent the distribution combined 
across all of the 5 fields; the dotted line is the distribution considering only the inner quarter of those fields, corresponding to distances out to about 700 kpc (at 
z ~ 1.5) from the AGN; the shaded region represents the expected field distribution, as derived at < 20 by combining the CADIS 16hr field (Thompson 
et al. 2000) and the K20 survey (Cimatti et al. 2002) and for 20 < if < 21 from the compilation of Hall & Green (1998). 



overdensities is not more sharply peaked should not be surprising. 
A galaxy cluster with a mass of a few x 10^* and an initial scale 
of 3-5 Mpc has an overall density of 1^ x 10^^^ kg m^^. The time 
that such a structure will take to separate out from Hubble flow and 
collapse is t ~ (37r/32Gp)^/^ ~ 3-6 Gyr (Peacock 1999). The 
Universe at redshifts z ~ 1.5 is about 4.5 Gyr old, and so such 
structures will still be collapsing into a virialised cluster such as 
those we see today; the large-scale structures seen around the AGN 
might even not yet be gravitationally bound. 

4.4 Bright red galaxies and sub-structure 

At high redshifts most clusters observed to date have several giant 
galaxies within them, presumably due to merging subclusters. Fur- 
ther, there is some evidence that AGN host galaxies are not always 
the only bright galaxies in their environment, for example due to 
sub-cluster mergers or being part of a group falling into a larger 
structure (e.g. Bremer, Baker & Lehnert 2002; Simpson & Rawl- 
ings 2002). Therefore, the fields have been searched for any galax- 
ies other than the AGN that are brighter than K = 18 and satisfy 
the colour criteria of the 'cluster elUptical candidates' above. The 
images around these bright red galaxies have then been examined; 



the majority of these galaxies appear to be fairly isolated, and so 
are likely to be foreground extremely red objects with redshifts of 
z ^ 1, since such foreground objects will tend to dominate the ex- 
tremely red galaxy population at brighter magnitudes (e.g. Cimatti 
et al. 2002). However, in two cases these bright red galaxies have 
associated fainter red galaxies, and may represent substructures at 
the AGN redshift. 

The first of these cases is in the field of 0000—177, where 
three galaxies with K = 17 Al, K = 17.77 wad K = 17.85 and 
colours — K > 5 are found close together, along with fainter 
red galaxies, about an arcminute east of the AGN. Indeed, these 
galaxies represent the small peak in the red galaxy radial distribu- 
tion (Figure 14) at 600 to 700 kpc radius around this source. This 
region of the field is of comparable richness in red galaxies to the 
regions around the AGN. The second case is a group found around 
?^ K = 17.09 galaxy about 1.5 arcminutes to the south-west of 
0016—129, partially causing the 800-1000 kpc bvunp in the radial 
distribution around that source. This appears even richer than the 
red galaxy distribution around the radio galaxy itself, which might 
explain why the value of iVo.s for this source is fairly low even 
though the number counts are well above the literature average. 

It is not clear whether these extra structures are indeed at the 
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Figure 14. The radial distribution of galaxies with 17.5 < K < 20.7 and R — K > Am each of the five fields separately, and the average of these. The radius 
is calculated in kpc at the redshift of each AGN. The error bars represent the Poissonian errors on the galaxy counts in each bin. In each figure the horizontal 
dotted line represents the expected blank-field distribution. 



same redshifts as the AGN, but if so then this would indicate that 
any large-scale structure around the AGN is still in the process of 
building up. It is interesting that this evidence for substructure is 
found in the two cases with the weakest evidence for red galaxy 
excesses in the 200-500 kpc radius range. 



low-concentration, irregular clusters. They interpret this as mean- 
ing that the mechanisms that produce morphological segregation 
(e.g. dynamical friction, or conversion of late-types into early types 
by gas stripping and halting star formation) work most rapidly in 
high mass systems. 



4.5 The morphology-density relation 

In nearby clusters, a well-known morphology-density relation ex- 
ists (Melnick & Sargent 1977; Dressier 1980), whereby in regions 
of high galaxy density, such as the centre of the cluster, higher pro- 
portions of early-type galaxies are found. Dressier et al. (1997) 
investigated the evolution of this morphology-density relation with 
redshift, out to a ^ 0.5, and find that this is still present at these 
redshifts but only in high-concentration, regular clusters; no corre- 
lation between galaxy-type and local galaxy density is found for 



The nature of the small-scale galaxy excess around these dis- 
tant radio sources (and in the two potential substructures) is par- 
ticularly interesting because, where present, it is comprised almost 
entirely of red galaxies (cf. Figs 4 to 6). This indicates that morpho- 
logical segregation may already have taken place in the very inner 
regions around the AGN, implying that the morphology-density 
relation is imprinted into cluster centres at a very early epoch. 
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4.6 Angular cross-correlation analyses 

The clustering of galaxies around the radio galaxies can be inves- 
tigated using the angular cross-correlation function w{9), which 
is defined from the probability (3P) of finding two sources in 
areas SQi and dfl2 separated by a distance 0: 5P = iV^[l + 
■w{9)]50,iSQ,2, where N is the mean surface density of sources on 
the sky. 

The value of w{9) can be evaluated for all galaxy-galaxy 
pairs, for just AGN-galaxy pairs, or restricting the galaxies in- 
cluded to certain ranges of colours and magnitudes. Here, four dif- 
ferent set of galaxy magnitudes and colours were considered: (1) 
'all galaxies'; (2) 'blue galaxies'; (3) 'red galaxies'; and (4) 'clus- 
ter elliptical candidates', which are the subset of the red galaxies 
that have K magnitudes, R — K and J — K colours in the range 
that would be expected for old passive cluster galaxies. The exact 
colour and magnitude definitions of each class of objects, and the 
number of galaxies contained in that class, are provided in Table 5. 

For each of these subsamples of galaxies the value of w{9) 
was estimated for a variety of bins in 9, following the method 
described by BOO, considering both all galaxy-galaxy pairs and 
only AGN-galaxy pairs. The results obtained over all 5 fields with 
multi-colour data were averaged to produce overall results, with 
an estimate of the uncertainty in the value of the combined w{9) 
in each bin provided by the scatter in the values between the 
different fields. w{d) usually shows a power-law form, w{9) = 
A{9/deg)~^, with a canonical value for the exponent of 6 = 0.8; 
this value appears also to be appropriate at high redshifts (Gi- 
avalisco et al. 1998), and so was adopted here for all fits here. This 
does not always provide an excellent fit to the data, but the current 
data do not have sufficient signal-to-noise to warrant an attempt to 
independently fit this parameter; note that because this parameter 
is fixed, the formal errors on the cross-correlation ampUtudes may 
be slightly higher than those quoted. 

The finite area of the images leads to a small correction, 
known as the integral constraint (C), such that w(6')obs ~ 
A[{9/ deg)^"'^ — C]. The value of C can be estimated by inte- 
grating w{6) over the area of each field, and for the current data 
corresponds to Cgg = 16.7 and Cag « 19.1 for all galaxy-galaxy 
pairs and just AGN-galaxy pairs respectively. Using these expres- 
sions, the amplitudes of the angular cross-correlation function were 
derived for both all galaxy-galaxy pairs and just AGN-galaxy pairs 
for each of the selected colour-magnitude subsamples of galaxies. 
The cross-correlation functions are displayed in Figure 15, and Ta- 
ble 5 provides the derived cross-correlation amplitudes. 

The angular cross-correlation amplitude of 'all galaxies' in 
these fields is within the range of the values recently derived to 
K = 20.75 for 'blank-fields' (e.g. Carlberg et al. 1997; Roche 
et al. 2002; Daddi et al. 2000a). However, comparing the values of 
Agg and /lag in Table 5 shows that the galaxies are much more 
highly clustered around the AGN than around other galaxies in 
general, indicating that a significant proportion of the galaxies are 
connected with the AGN. Evidence for such clustering around the 
AGN is particularly apparent when comparing the results for galax- 
ies of different colours. The blue galaxies show no evidence for 
clustering, not even around the AGN, whilst red galaxies are known 
to be more highly clustered (e.g. Roche et al. 1996, Shepherd et al. 
2001 and references therein) and this is reflected in the value of 
Agg derived for this class; their higher value of .4ag shows that 
these are preferentially clustered around the AGN. The 'cluster el- 
liptical candidates' have even larger values of Agg and Aa,g, sug- 
gesting that this selection criteria picks out a highly clustered class 



of object, many of which are indeed likely to be old passive cluster 
ellipticals. It should be noted that this clustering signal is domi- 
nated by the points at small radii, < 10 arcsec, corresponding to 
the excess in the upper left panel of Figure 12; however, even if 
the inner two bins of the cross-correlation analysis are ignored, the 
clustering amplitude around the AGN is still higher than that of the 
field in general, indicating that it is not only these nearby galaxies 
which produce the clustering signal. 

If a population of A^c cluster galaxies, with a high intrinsic an- 
gular cross-correlation amplitude (Ac), is observed together with 
a sample of M foreground or background interloper galaxies with 
zero clustering amplitude, then the resultant combined sample will 
have an observed cross-correlation amplitude of ~ j^^^^Ac. In 
other words, the observed amplitude will be the amplitude that 
would be found for the cluster galaxies alone, scaled down by the 
proportion of the total sample that are cluster galaxies. This sit- 
uation is a good approximation to those of the 'red galaxy' and 
'cluster elliptical candidate' samples, since the observed value of 
Agg for these samples are almost an order of magnitude below 
those of Aag. Comparing the values of ^ag then indicates that the 
colour definition for the 'cluster elliptical candidate' sample ap- 
proximately halves the fraction of interlopers relative to the red 
galaxies. 

If the results from the number count analysis are taken at face 
value, then there are about 270 clusters galaxies across the five 
fields studied in the cross-correlation analysis (290 excess galax- 
ies determined in Section 4.1, of which 20 are from the 2128—208 
field, not considered here); using this, an estimate can be made for 
the number of cluster galaxies in each of the colour-selected sam- 
ples. Assuming that about 80% of the cluster galaxies are 'red' 
(from Figure 13) then within the 'red galaxy' sample there are 
about 220 cluster galaxies, which corresponds to about 25% of the 
overall 'red galaxy' sample (see Table 5). This would suggest that 
about 40% of the 'cluster elliptical candidate' sample, that is ~100 
galaxies, are indeed cluster galaxies. If these results are correct, the 
~100 early-type cluster galaxies at redshifts z ^ 1.6 would pro- 
vide a substantial sample for cluster galaxy evolution studies. 

As has been described by many authors (e.g. Longair & Seld- 
ner 1979, Prestage and Peacock 1988), if the form of the galaxy lu- 
minosity function at the redshift of the AGN is assumed, it is possi- 
ble to convert the amplitude of the angular cross-correlation func- 
tion for galaxies around the AGN into a spatial cross-correlation 
amplitude, Bag, derived from the spatial cross-correlation function 
^(r) = Bag ( ■ Following the method of BOO, the value of 
Aag = 0.0061 ± 0.0022 for the 'all galaxies' sample corresponds 
to a spatial cross-correlation amplitude of Bag = 660 ± 240. This 
amplitude is comparable to the value of Bag = 510 ± 120 derived 
by BOO around 3CR radio galaxies with redshifts z ~ 1. These 
values can be interpreted physically by comparing with the equiva- 
lent values (Beg) for Abell clusters calculated between the central 
galaxy and the surroimding galaxies; these have been derived by 
several authors and average to Beg « 350 for Abell class and 
Beg « 710 for Abell class 1 (Prestage & Peacock 1988; Hill & 
Lilly 1991; Andersen & Owen 1994; Yee & Lopez-Cruz 1999). 
The environments surrounding these z ~ 1.6 AGN seem to be, on 
average, comparable to those of local clusters of richness between 
Abell classes and 1. 
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Figure 15. The angular cross-correlation functions for galaxies with a range of magnitudes and colours as defined in Table 5, averaged across the 5 fields. The 
left-hand plots correspond to all galaxy-galaxy pairs and the right-hand plots consider only AGN-galaxy pairs. 



Table 5. Results from the angular cross correlation analyses. For galaxies with a range of magnitudes and colours, the amplitude of the angular cross-correlation 
fimction is provided considering both all galaxy-galaxy pairs (Agg) and only AGN-galaxy pairs (Aag). 

Angular cross-correlation amphtude 



Sample 


K-mag range 


IR col range 


IR-opt col range 




All galaxies 


K<20.7 


All 


All 


1739 


Blue galaxies 


K<20.7 


J-K<1.6 


R-K2.5 


762 


Red galaxies 


K<20.7 


J-K>1.4 


R-J>2.0 


868 


Cluster ellips? 


17.0<K<20.5 


1.6<J-K<2.2 


3.0<R-J<4.5 


252 



0.0005 ± 0.0001 
0.0000 ±0.0003 
0.0020 ± 0.0002 
0.0042 ± 0.0004 



'lag 

0.0065 ± 0.0022 
0.0000 ± 0.0026 
0.0166 ±0.0041 
0.0307 ± 0.0052 



5 ENVIRONMENTAL VARIATIONS COMPARED WITH 
RADIO SOURCE PROPERTIES 

The results have clearly shown that the environments of these high 
redshift AGN are heterogeneous, with considerable variation in 
both the amplitude and the physical extent of any galaxy excess. 
In this section the variations in the galaxy overdensities are com- 
pared with the other properties of the radio source, to see if there 
are any indications of why the enviroimients differ so much. The 
properties of the radio sources studied are provided in Table 6. 

The first thing to note is that all of these radio sources have 
essentially the same redshifts and radio powers, and so the varia- 
tions are not generally correlated with either of these properties. 
Equally, all 6 studied sources are steep-spectrum radio sources ac- 
cording to the standard definition of a > 0.5 (where a is defined as 
Sv oc v~"), although 1524—136 has a spectral index close to this 
cut-off; the weaker small-scale excess of this source could con- 
ceivable be related to this in some way. No fundamental difference 
is seen between the quasar and radio galaxy populations, and nor is 



any correlation with the size of the radio source seen for any of the 
extended radio sources. 

The only issue which this comparison does bring to light is 
that the one source for which no overdensity is seen on either 
large or small scales, 2128—208, is an unresolved radio source; this 
could potentially explain this result. Extended radio sources must 
have been fueled continuously for timescales of at least \(f years, 
and possibly up to 10^ years, which requires a plentiful and steady 
supply of fueling gas, as well as a relatively dense surrounding 
medium to confine the expanding radio lobes. Perhaps at high red- 
shifts these conditions are only satisfied for galaxies within larger- 
scale structures. On the other hand, compact radio sources are rel- 
atively young, < 10 ' years old (e.g. Owsianik and Conway 1998), 
and it has been suggested that some proportion of these may fiz- 
zle out before evolving into larger radio sources (e.g. Alexander 
2000). The reduced gas supply required for these short-lived com- 
pact sources may then be available to galaxies in any enviroimient, 
and could explain the lack of any overdensities arotmd 2128—208. 

This situation has been modelled by Kauffmann & Haehnelt 
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(2002) who have investigated the quasar-galaxy cross-correlation 
function around quasars at different redshifts within hierarchical 
clustering models of galaxy formation. In their model they assume 
that quasars are triggered by galaxy mergers, with some percentage 
of the gas brought in by the merger then fuelling the black hole. 
They show that the cross-correlation function around the quasars is 
significantly higher for quasars with long lifetimes than for short- 
lived sources, in agreement with the discussion above. This differ- 
ence is most pronounced for the most powerful AGN, such as those 
studied here. 

An alternative may be that this is a lower power radio source 
with its fltix density Doppler boosted by beaming; in this case a 
correlation between radio power and environmental density would 
explain the result. However, a beaming scenario is considered un- 
likely due to the steep spectral index of the source. 



6 DISCUSSION AND CONCLUSIONS 

The results of this work can be summarised as follows. 

• On average, a significant overdensity of K-band selected 
galaxies with 17.5 < K < 20.5 is found across the AGN fields. 

• These excess counts are predominantly associated with galax- 
ies with magnitudes K > 17.5 and colours R — K > A, which 
are (although not uniquely) exactly those expected for old passive 
cluster galaxies at redshifts « ^ 1.6. 

• These galaxies are fovmd on two different spatial scales arovmd 
the AGN: pronounced concentrations of galaxies at radial distances 
< 150 kpc, and weaker large-scale overdensities extending out to 
between 1 and 1.5 Mpc. 

• The presence or absence of galaxy excesses on these two 
scales differs greatly from AGN field to AGN field. All of the ex- 
tended radio sources show overdensities on some scale: two show 
overdensities on both scales, two predominantly on large scales and 
one only on small scales. The one field which shows little evidence 
for a galaxy excess on any scale is associated with an unresolved 
radio source. 

• Where overdensities are present on < 150 kpc scales, these 
are composed almost entirely of red galaxies. This suggests that 
the morphology-density relation is imprinted into the centres of 
clusters at a very early stage of cluster formation. 

• The amplitude of the angular cross-correlation function 
around the AGN is a strong function of galaxy colour. The selection 
of only galaxies with magnitudes and colours consistent with being 
old passive elliptical galaxies at the AGN redshifts maximises this 
amplitude, indicating that old ellipticals do indeed exist in these 
environments. 

• The spatial cross-correlation amplitude for all galaxies around 
the AGN is similar to that around the 3CR radio galaxies at red- 
shifts z ^ 1; both are consistent with having the same richness as 
local Abell clusters of richness class to 1. 

It is interesting to consider what these results imply for the nature 
of the AGN at these redshifts, and also for cluster evolution. 

Powerful radio galaxies are known to be hosted by giant el- 
liptical galaxies at all redshifts z < 1.3 (e.g. Best et al. 1998, 
McLure & Dunlop 2000), and even out to redshifts « < 2.4 a 
significant proportion still have radial profiles well-matched by de 
Vaucouleurs' law (Pentericci et al. 2001). Little difference is found 
between the host galaxies of radio galaxies and radio-loud quasars 
at these redshifts (Kukula et al. 2001), which suggests that the hosts 



of all of the 2 ~ 1.5 AGN studied here are likely to be giant ellip- 
tical galaxies. 

The variety of different environments fotmd for these AGN 
may simply reflect this fact: that the primary requirement for form- 
ing a powerful radio-loud AGN is a giant elliptical galaxy host. 
Such galaxies can be found in a range of environments from galaxy 
clusters, through groups of galaxies, to relatively isolated environ- 
ments, but are generally only formed by redshifts 2; ~ 1.6 when 
they are found at the highest peaks of the primordial density per- 
turbations, within larger-scale structures. This explains the high 
average richness, but still significant variation, of the large-scale 
galactic environments arovmd these AGN. 

To understand the pronounced smaller scale galaxy overden- 
sity, the key question is whether these central galaxies represent 
a stable virialised structure, or whether they are still infalling and 
will merge with the central radio galaxy. If these are already viri- 
alised then the small scale galaxy overdensity is clearly related to 
the larger-scale structure, and indicates that at least some of the 
AGN lie at a special location in their environments: at the centres of 
their clusters or groups. This again is not too surprising: the kinetic 
energy of the relativistic radio jets of these AGN corresponds to the 
Eddington luminosity of a black hole with M ~ 10^ M© (Rawl- 
ings & Saunders 1991), and so even assuming that these sources are 
fueled close to the Eddington limit they must be powered by very 
massive central engines and the radio lobes must be confined by 
dense gas. It is not unreasonable to expect these most powerful ra- 
dio sources to be found within the most massive elliptical galaxies, 
at the centre of the cluster or group of galaxies. 

In this interpretation, the four different scenarios of the pres- 
ence or absence of galaxy overdensities on small and large scales 
correspond to four different possible locations of the AGN host 
galaxy: it may be a central cluster galaxy in cluster with a viri- 
alised core (excesses on both scales), an off-centre giant elliptical 
in a cluster, or the central galaxy in a cluster at an early stage of 
formation (only large-scale overdensity), the central galaxy of a 
compact group (only small-scale excess), or an isolated giant el- 
liptical (no excess on any scale). With respect to this model, it is 
interesting that all of the extended radio sources lie in group or 
cluster environments. 

Alternatively, if the galaxies producing the small-scale excess 
are not yet virialised, then they may instead be related to a hier- 
archical build up of the central AGN host galaxy. Pentericci et al. 
(2001) found that the optical characteristic sizes of powerful radio 
galaxies are typically a factor of two smaller at 2; ~ 2 than those 
at « ~ 1. Since the galaxy luminosity and the characteristic size of 
ellipticals are related by Lint oc r|?'^ (Kormendy 1977), this means 
that the radio galaxies must grow in mass by about a factor of two 
over this redshift interval, a large part of which may be due to merg- 
ers with companion galaxies. If the 150 kpc scale galaxy excess has 
a velocity dispersion a ~ 500 km s^^, then the galaxy crossing 
time for this structure is ~ 0.3 Gyr. A galaxy merger timescale of a 
few crossing times is therefore comparable to the interval of cosmic 
time between redshifts 2; ~ 2 and 2 ~ 1, and some of the small- 
scale excesses of galaxies may be destined to merge with the AGN 
host galaxy soon. 

Whichever of these interpretations for the small-scale galaxy 
excess aroimd the AGN is correct (and undoubtedly a combina- 
tion of both effects is occurring), the Mpc-scale galaxy overdensi- 
ties appear be associated with (proto?) cluster environments. The 
wealth of galaxies with magnitudes and colours consistent with be- 
ing old elUptical galaxies at 2 ~ 1.6 indicate that the epoch of 
elliptical galaxy formation must be earlier than this redshift; the 
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Table 6. Properties of the 6 radio sources studied. 



Source 


z 


*S'408 


„408 
"1400 






[Jy] 




MRC 0000-177 


1.47 


6.51 


0.80 


MRC 0016-129 


1.59 


6.87 


0.95 


MRC 0139-273 


1.44 


5.04 


0.95 


MRC 1524-136 


1.69 


6.11 


0.61 


MRC 2025-155 


1.50 


5.41 


1.05 


MRC 2128-208 


1.62 


6.15 


0.88 



* 'Y' = strong-excess; 'y' = 
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^^408 


D 


Type 


Excess on scale* 


[lO^^'WHz-i] 


["] 




<150kpc 




85.4 


2.7 


Q 


Y 


N 


125.4 


3.5 


RG 


Y 


y 


72.2 


12 


RG 


y 


Y 


92.5 


0.4 


Q 


y? 


Y 


93.9 


15 


Q 


Y 


y 


110.1 


< 1 


Q 


N 


N 


weak excess; N = 


no significant excess. 





colours of these objects further suggest that the bulk of their stellar 
populations must have formed at redshifts 2 > 5 (this will be dis- 
cussed further in Best et al. in preparation). By using these radio- 
loud AGN as signposts, it is therefore possible to define significant 
samples of z ^ 1.6 cluster early-type galaxies with which to ex- 
tend studies of galaxy evolution back to earlier cosmic epochs. 

Finally, it is interesting to compare the comoving number den- 
sity of these powerful radio sources with that of nearby rich clus- 
ters. According to the pure luminosity evolution model for steep- 
spectrum radio sources of Dunlop and Peacock (1990), at redshifts 
z = 1 to 2, there are about 10~^ radio sources per comoving 
Mpc^ within the upper decade of radio luminosities containing the 
sources studied here. By comparison, the comoving number density 
of nearby clusters of Abell richness class or greater is ~ 4 x 10~® 
Mpc^'\ and for Abcll class 1 and richer it is ^ 1.7 x 10^'' Mpc"^ 
(e.g. Bahcall & Soniera 1983; converted to the cosmology adopted 
in this paper; cf. West 1994). However, radio sources only live for 
between lO'^ and 10* years (e.g. Kaiser et al. 1997), and so the 
observed number density of powerful radio sources underestimates 
the number density of environments which will contain a powerful 
radio source at some point between redshifts 2 and 1 by a factor 
~ 50 (= cosmic time between redshift 2 and 1 divided by the life- 
time of a typical radio source). Therefore, the comoving number 
density of environments which will contain a powerful radio source 
between redshifts 2 and 1 is comparable to the comoving number 
density of nearby Abell clusters of richness class or greater. This 
may simply be an interesting coincidence, or alternatively may im- 
ply that the onset of a powerful radio source is a key stage in the 
formation of a rich cluster. 
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